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Observation of D — tt Production Correlations 
in 500 GeV tt^ — N Interactions 
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Abstract 



We study the charge correlations between charm mesons produced 
in 500 GeV tt" — N interactions and the charged pions produced closest 
to them in phase space. With 110,000 fully reconstructed D mesons 
from experiment E791 at Fermilab, the correlations are studied as func- 
tions of the Dtt — D mass difference and of Feynman x. We observe 
significant correlations which appear to originate from a combination of 
sources including fragmentation dynamics, resonant decays, and charge 
of the beam. 
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While the production of heavy quarks can be calculated in perturbative Quantum 
Chromodynamics (QCD), the evolution of these quarks into hadrons remains one of 
the most challenging aspects of nonperturbative QCD. Correlations between charm 
mesons and the charged pions produced closest to them in phase space provide infor- 
mation on how quarks evolve into hadrons. Fragmentation dynamics 0], resonances 
0, and beam effects can each produce such correlations. The relative importance of 
these mechanisms must be determined from data. 

During fragmentation, correlations could be produced because qq pairs from the 
vacuum are neutral. For example, if a c quark combines with a d from such a pair to 
form a , the remaining d is close by in phase space and is likely to become part 
of the closest pion, which we call the "associated pion". Thus, D^n^ {D^n^) would 
be favored and D^tt^ {D^tt^) disfavored. Similarly, D^n^ {D^n^) would be favored 
and D^TT^ (D^-k^) disfavored. Resonances produce the same favored associations. 
D*^ decay associates a tt^ with a while D*~ decay associates a vr" with a D^. 
Qualitatively, D** decays produce the same correlations. 

The charge of the beam particle can also lead to charge correlations. Using a vr^ 
beam can lead to the association of both charm mesons and anticharm mesons with 
negative pions, especially in the forward (beam) direction. Two distinct but related 
mechanisms can lead to this result. If the charm quark (antiquark) produced in a 
hard interaction coalesces with the antiquark (quark) from the beam particle to form 
the charm (anticharm) meson, the remaining quark (antiquark) from the beam can 
become part of a negative pion, but not part of a positive pion. If neither the quark 
nor the antiquark from the beam pion is used in making the charm meson, both are 
available to form negative pions but not positive pions. 

By comparing the charge correlations of different species of charm mesons and 
antimesons with associated pions, and by studying them as functions of Feynman 
X (xp) of the charm meson, one can hope to disentangle some of these processes. 
Evidence of such correlations between B mesons and associated light mesons, ascribed 
to resonances, has been observed in decays at LEP by the OPAL collaboration 
1^. In this letter, we report the first observation of fragmentation and beam- related 
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production correlations for charm mesons. 

We use D° — i> _ft'^7r+, K^n^n^ , and D*^ D°7r+ signals (and their charge 

conjugate decays) from experiment E791 at Fermilab for this study. The data were 
recorded using a 500 GeV/c vr^ beam interacting in five thin target foils (one platinum, 
four diamond) separated by gaps of about 1.4 cm. The detector, described elsewhere 
in more detail is a large-acceptance, forward, two-magnet spectrometer. Its key 
components for this study include eight planes of multiwire proportional chambers, 
six planes of silicon microstrip detectors (SMD) before the target for beam tracking, 
a 17-plane SMD system and 35 drift chamber planes downstream of the target for 
track and vertex reconstruction, and two multicell threshold Cerenkov counters for 
charged particle identification. 

During event reconstruction, all events with evidence of well-separated production 
and decay vertices were retained as charm decay candidates. For this study, we require 
the candidate charm decay vertex to be located well outside the target foils and to be 
at least 8(Ta downstream of the primary vertex (where cta is the error in the measured 
longitudinal separation between the vertices ~ 350/im). The momentum vector of 
the candidate D must point back to the primary vertex with impact parameter less 
than 80/im. The momentum of the D transverse to the line joining the primary and 
secondary vertices must be less than 0.35 GeV/c. Each decay track must pass closer to 
the secondary vertex than to the primary vertex. Finally, the track assigned to be the 
kaon in the charm decay must have a signature in the Cerenkov counters consistent 
with the kaon hypothesis. The D*^ candidates are found from the / samples 
by adding vr^ tracks and requiring that the mass difference Am = M{Dtx) — M{D) 
be consistent with the D* Dn hypothesis. The final signal sizes are obtained by 
fitting the invariant mass spectra as Gaussian signals and linear backgrounds. For 
W, D+, D~, D*+, and D*" , the fits yield 22587 ±210, 24237 ±216, 24569 ±204, 
29649±238, 4997±84 and 6048±93 events, respectively. The r.m.s. mass resolutions, 
aD, used later in defining signal and background bands, are 13 MeV/c^, 13 MeV/c^, 
and 14 MeV/c^ for D^, and respectively. 

For each D found in an event, all tracks originating from the primary vertex 
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and producing a pion signature in the Cerenkov counters are combined with the D. 
Among these combinations, the pion that forms the smallest invariant mass difference 
(Ammm) with the D decay products is selected as the associated pion. 

We define the correlation parameter a as 

where g = +, -, -, +, -, + for D = , D° , D+, D~, D*^ , and D*~ , respectively, 
and Ni{DTT'^) denotes the number of charm mesons for which the selected pion has 
the charge q. In the absence of correlations a is zero, and in maximally correlated 
cases it is unity. 

We first study the Dn correlations as functions of Ammin for Arrimin < 0.74 
GeV/c^. The number of Dtc signal combinations in each Arrimin bin is determined 
by subtracting from the Arrimin distribution for D candidates (mass within ±2.5 ctd 
of the nominal D mass) the appropriately normalized Arrimin distribution for back- 
ground events (mass between 3.0 ctd and 5.5 from the nominal D mass). The 
correlation parameters for background-subtracted signals (before additional correc- 
tions) and background regions are listed in Table 1. The signal correlations differ 
significantly from the background correlations. We note that replacing the D can- 
didate in an event with a D of the same species from another event, while keeping 
the rest of the event the same, produces correlations consistent with those of the 
background. 

We use a Monte Carlo simulation of the experiment and the LUND event generator 
(PYTHIA 5.7/JETSET 7.3) [|] to model the effects of our apparatus and reconstruc- 
tion. This simulation describes the geometry, resolution, noise, and efficiency of all 
detectors, as well as interactions and decays in the spectrometer. The detected D*/D 
production ratio in the Monte Carlo matches our data well. As with real events, the 
associated pion for each reconstructed D meson is selected. By matching the selected 
pion's momentum vector with the momenta of all generated particles, we determine 
whether the selected pion track is a real track or a ghost (false) track. Selecting a 
ghost pion (not matched to any generated track) or a real pion not matched to the 
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true associated pion can cause smearing in Arrimin and dilution of the correlation. 
Selecting a pion with the same charge as the associated pion but with different mo- 
mentum smears events in Am^j„. Selecting a pion with the opposite charge smears 
events in Am^in and also dilutes the correlation. 

To account for the effects of ghost tracks, smearing, dilution, and acceptance on 
the correlations as functions of Am^j„, we employ a matrix formalism. For the 
the observed number of D'^tt^ combinations 0~j^ in the bin of Am^jn can be 
written as 

or = E S]fArNt^ + S^fAl^Nt^ + G^Oj^ (2) 

i i 

where N^^ denotes the true number of D^-k^ events in the ith bin of Arrimin, Af^ 
the acceptance probability, and G'j'^ the ghost track rate for D'^tt'^ combinations. 

The matrix S^^ describes smearing in the absence of dilution while the matrix S'^^ 

describes smearing and dilution when the wrong sign pion is selected. The smearing 

matrices S^^ and 5"^^, the acceptance coefficients and A'^'^, and the ghost track 

rates G'^~ and are determined from the Monte Carlo. The coupled matrix 

equations in (2) are solved to obtain the true distributions and N^'^. Corrected 

Ammm distributions are shown in Figure 1. The corrected correlation parameters for 

D, a{D), for D, a{D), and for the D and D combined, a{D,D) are presented in 

column 4 of Table 1. 

The statistical and systematic errors assigned to the final measurements, shown 
first and second respectively, are also given in Table 1. These errors are propagated 
through the matrix formalism. The systematic errors account for uncertainties in 
the Monte Carlo simulation of the detector (their effects on dilution, smearing, ghost 
tracks, and acceptance), analysis cuts, background subtraction, kaon misidentifica- 
tion, and binning (in decreasing order of importance as listed). For each data point, 
the systematic uncertainties due to these sources are added in quadrature. The sys- 
tematic uncertainties due to statistical fluctuations in the Monte Carlo are neghgible. 

To verify the results produced by the matrix formalism, we also estimate the 
correlations using simple dilution factors (summed over all bins of Arrimin)- For D~^, 
the true number of combinations, N^~~ and N^"^, can be expressed in terms of the 
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reconstructed combinations iV+ and N^~^ as 

iV+^ = (1 - rf+^)iV+^ + rf+±iV+±, (3) 

where the dilution factor d^^ denotes the probabihty that a true D^n^ combination 
is reconstructed as a D+vr^. The results from this technique are consistent with those 
reported in Table 1. 

All studies and corrections have been done within the framework of the LUND 
PYTHIA/JETSET model. The dilution factors dab in Eq- (3) are typically of order 

0.2 — 0.3. In our Monte Carlo, ~ but d_^ is less than d The difference 

between d and d ^. is almost independent of xp with a typical value near 0.06. 

Varying some of the JETSET fragmentation parameters to reproduce our inclusive 

/ production asymmetries as a function of xp, as described in ref. 0, leads to 
results consistent with those in Table 1. A fundamentally different model of hadron 
production might change the differences between the d's discussed above by a few 
times 0.01, which would in turn change the measured correlation parameters. For 

example, reducing {d — d ^) from 0.06 to 0.05 would increase a{D^) by 0.02 - 

0.03. 

In Figs. 1(a) and (b) we present the numbers of D^tt^ and D^ti^ combinations 
as functions of Arrimin- In both of these plots the combinations differ mainly in the 
D*^ resonance region ,the first 75MeV/c^ bin. Using a ±2.5(T cut on the D*^ — 
and D*^ — mass difference, we separate the final D^tc^ and D^n~ samples into 
resonance (res) and continuum (cont) contributions to obtain a{D^^g) = 0.98 ± 0.04 
and a{D^res) = 0.98 ± 0.02. For pure resonance, a would be near 1. The measured 
values serve as a check of our method. The continuum measurements are a^D^^^^) = 
-0.07 ± 0.03 and a{W^ont) = 0.17 ± 0.03. In Figs. 1(c) and (d) we present the 
D^Ti^ and D^n^ combinations. In both these plots the combinations differ over a 
broad range in Am^m- In Figs. 1(e) and (f) we present the D*^7i^ and D*^7T^ 
combinations. A pattern similar to that for is manifest. The plots for charm 
mesons and anticharm mesons clearly differ. These differences also appear in column 
(4) of Table 1, and indicate the presence of significant beam-related effects. 

To investigate beam-related effects in more detail, we study the xp dependence of 
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the and D*^ correlations. We do not show the correlations for I?°'s since many 
£)°'s are decay products of either D*^ and making interpretation difficult. In 
Fig. 2, we plot q; as a function of xp, for both particle and antiparticle for and 
D*^ . The distributions are corrected using the simple dilution factor technique. We 
observe that a{D^) rises shghtly with xp but Oi{D~) falls sharply to negative values 
for xp > 0.2. In both cases, the D^s are more hkely to be associated with 7r~'s at high 
xp where beam effects seem to be important. There appears to be less dependence 
of a on a^i;' for the D*"^. 

Further beam-related studies with Monte Carlo data suggest the correlation asym- 
metries cancel under neutral beam conditions and are in fact symmetric. This is 
effectively accomplished when the combined particle and antiparticle correlations are 
computed. In Table 1, wc show the combined and symmetrized correlation param- 
eters to be a{D^,W) = 0.29 ± 0.02 ± 0.03, a{D+,D-) = 0.21 ± 0.02 ± 0.03, and 
D*~) = 0.23 ± 0.04 ± 0.03. These results indicate that fragmentation dynam- 
ics and resonant decays produce substantial correlations between D mesons and their 
associated pions. All three combined correlation levels are approximately equal, al- 
though the correlations for neutral and charged D mesons are dominated by resonant 
and continuum regions of Arrimin, respectively. 

In summary, we observe significant production correlations between D mesons and 
their associated pions. Some of these correlations are associated with fragmentation 
dynamics, some with resonances, and some with the charge of the beam. In addition 
to providing information on how heavy quarks evolve into hadrons, such correlations 
may provide tools for tagging flavor in CP violation studies in heavy flavor systems. 
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TABLES 

TABLE I. The xp- and Am^i„-integrated correlation parameters a defined in Eq. (1) 
for the background-subtracted signals prior to correction, for the corresponding back- 
grounds, and for the signals after correction using the matrix technique based on Eq. (2). 
a{D^) and a{D^) contain both D* resonance and nonresonance contributions. 



Charm 


Signal a 


Background a 


Corrected Signal a 




0.13 ±0.01 


-0.04 ±0.01 


0.12 ±0.03 ±0.04 




0.18 ±0.01 


0.04 ± 0.01 


0.42 ± 0.02 ± 0.03 




0.16 ±0.01 


0.00 ± 0.01 


0.29 ± 0.02 ± 0.03 


D+ 


0.18 ±0.01 


0.10 ±0.01 


0.45 ± 0.03 ± 0.03 


D- 


0.08 ± 0.01 


0.02 ± 0.01 


0.03 ± 0.03 ± 0.04 


D+,D- 


0.13 ±0.01 


0.05 ±0.01 


0.21 ±0.02 ±0.03 


D*+ 


0.15 ±0.02 


0.08 ± 0.03 


0.33 ± 0.05 ± 0.03 




0.08 ± 0.02 


0.02 ± 0.03 


0.15 ±0.05 ±0.04 


D*+,D*- 


0.12 ±0.01 


0.05 ± 0.02 


0.23 ± 0.04 ± 0.03 
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FIGURES 
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FIG. 1. The fully corrected Arrimin distributions for {a)D^n^, {h)D^Tr^, (c)L>+7r'F, 
{d)D-Tr^, (e)L'*+7r'F, and (f)i:>*-7r=^ combinations. The error bars are statistical only. 
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FIG. 2. The corrected correlation parameter a as a function of xp for {a)D~^ and 
(h)D*^ . The parameter a is defined in Eq. (1) in the text. Tlic error bars correspond to the 
statistical and systematic uncertainties added in quadrature. Additional model-dependence 
is discussed in the text. 
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